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Progress report NAGW-3954 

The objective of this program is to quantify, by laboratory experiments, the charging of ices and 
other insulators subject to irradiation with electrons, ions and ultraviolet photons and to model 
special conditions based on the data. The system and conditions to be studied are those relevant 
for charging of dust in magnetospheric plasmas. The measurements are supplemented by 
computer simulations of charging of grains under a variety of conditions. 

Our work for this period involved experiments on water ice, improved models of charging of ice 
grains in Saturn's E-ring, and the construction of apparatus for electron impact studies and 
measurements of electron energy distributions. 

Electron emission and surface potentials of ice under ion bombardment 

We studied electron emission and sputtering during irradiation of amorphous water ice at 60 K by 
He*, Li*, Be* B*, C*, N*, O* F*, and Ne* ions in the energy range from 10 to 100 keV. 

The electron yields where determined with the standard charge-collection method using suitably 
biased electrodes [1,2] and the sputtering yields with an ultra sensitive quartz-crystal 
microbalance that can detect the removal of a tenth of a monolayers of water [3], 

We found that for constant velocity (5 keV/amu) ions the dependence of the sputtering yields 
with projectile atomic number Z is proportional to the square of the Z-dependence of the 
electronic stopping power dE/dx. The electron yields increase sublinearly with dE/dx, contrary to 
assumptions made so far in the literature, which just applied results for metal targets [4,5], We 
explain this sublinear effect by noting that unlike in metals, where electrons are very mobile and 
can neutralize any excess charge, ionizations produced in insulators like water ice cannot be 
screened during electron emission, due to the localization of the electrons. This means that the 
ionization track produced by each projectile will charge up, transiently, to a potential that will act 
as a barrier for electron escape. This is a localized charging, that superimposes on any 
macroscopic charging [6], Contrary to expectations, this charging of the track is important for 
the weakly ionizing protons meaning that it will also be important for energetic electrons. 

The sputtering studies have shown that the yields are higher than what has been assumed in the 
past. This has important consequences on estimates of the source of plasma in the inner Saturnian 
magnetosphere [7], 

In another study, supported by NSF but relevant to this work, we found what we term 
"catastrophic sputtering" for condensed Ar: an increase of the sputtering yield, by up to two 
orders of magnitude, due to electrical discharges or breakdown produced by irradiation in a weak 
electric field. So far, we have not found this effect in pure water ice under normal conditions but 
we cannot rule it out for films thicker than the ion range, for other insulators, or for high 
environmental electric fields. 
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New experimental setup for electron irradiations 

We obtained a Perkin-Elmer Auger spectrometer which we are adapting to this work, having built 
already the alignment and support hardware and a differential pumping chamber to connect it to 
our 60 kV ion accelerator. With this system, we will be able to irradiate insulators with low 
energy electrons and keV ions and measure the energy distribution of the ejected electrons. 

Modeling the charging of ice grains * 

This year, we modified our previous modeling of charging of micron sized ice grains, for 
application to Saturn's E-ring. We took into account two aspects ignored in previous studies, 
which we found to be essential for grains in low temperature plasmas: electron reflection and the 
threshold in secondary electron emission due to the band gap of the insulator. We calculated the 
electrical potential of the grain from the balance between different fluxes of charged particles, 
their energy distribution, and the secondary emission properties of ice. Since the energy 
distributions and fluxes of incoming and outgoing particles are affected by the potential of the 
grain, the problem of dynamical charging is a complex one that needed to be treated self- 
consist ently. We found that existing charging models cannot be simply extrapolated to the low 
energies (<30 eV) common in planetary magnetospheres. Using parameters from the Voyager 
PLS experiment, we calculated the potential of Saturn's E-ring grains to vary from -5.5V to +5V 
at distances from 4 to 10 Saturn radii. 

Work for the second year 

The plan for the next year is to: 1 ) extend the results obtained using ions on water ice to other 
ices (CO2, CH4, NH3, etc.), 2) finish installation of the electron irradiation chamber and do the 
first experiments on secondary electron emission and energy distributions; 3) determine 
composition changes induced by bombardment and how they affect electron emission and 
sputtering properties of the ices. 

We will finish the installation of the new chamber for electron and ion bombardment. This 
involves: a) Constructing and installing cooled target stage with temperature control in the range 
20 K - 200 K; b) A gas manifold and gas-doser using a micro-capillary array, to grow uniform ice 
layers; c) Installing electronics and adapting the required software for data acquisition, d) building 
and installing electronics for use with pulsed projectile beams, to measure both the electron yields 
and the charging of the target. 

The first experiments will measure secondary electron yields. We will study the dependence on 
film thickness, amd focus on trying to achieve electron energies below 30 eV, which our modeling 
has show to be those which determine the grain potential in the E-ring environment. 

The measurements will also be done as a function of irradiation dose, since insulators typically 
decompose under impact of ionizing radiation [8], This decomposition may be accompanied by 
the loss of different components, at rates that depend on target temperature and which we can 
monitor by mass spectrometry, as shown in our recent work on ice using Lyman-a radiation [9]. 
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In other materials like methane, irradiation leads to the formation of less volatile compounds, as 
we have recently verified in our laboratory. To understand these processes, we will look at 
irradiation induced changes in the solid using photoelectron spectroscopy, by adding cooling 
capabilities to our Perldn-Elmer 560 ESCA surface analysis system. We will correlate in-situ 
measurements of composition changes induced by bombardment, and changes in electron 
emission properties. 

Computer simulations 

Our modeling has shown to be very valuable in establishing the conditions most important for 
charging and thus guide fixture experiments. We will start simulating the effect of grain size and 
topography by using Monte-Carlo simulations of electron emission in different parts of the grain 
which results when it is traversed by a fast projectile [10], This work will form the foundation of 
subsequent simulations that will address the complexities of varying surface compositions within a 
given grain, different grain shapes, thickness of ice mantles, etc. 

List of papers in the last year 
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Abstract 

The charging of ice grains in planetary plasmas is studied, including the 
effects of secondary electron emission and backscattering of the incident elec- 
trons. It is shown that existing charging models can not be simply extrapolated 
to the low energy electron regime (below 30 eV) common in planetary magne- 
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tospheric plasmas. We derive expressions for the electrical potential of a grain 
immersed in a low energy plasma which more carefully account for electron 
reflection and the threshold for secondary electron emission. Using plasma pa- 
rameters from Voyager PLS experiment, we calculate the potential of Saturn’s 
E-ring grains to vary from -5.5 V at 4 R, to 5 V at 10 R 3 . 

1. INTRODUCTION 

The E-ring is a diffuse, azimuthally symmetric distribution of small water-ice grains 
in Saturn’s magnetosphere occupying the region between 3 and 8 Saturnian radii. 
This ring appears to be composed predominantly of ice grains (1 ± 0.3 micrometer 
in radius, Showalter et al., 1991), which are surrounded by a low-density plasma 
consisting of electrons, protons and single-ionized oxygen ions ( Richardson and Sittler, 
1990). Three spacecraft have traversed Saturn’s magnetosphere and measured plasma 
parameters: Pioneer 11, Voyager 1 and Voyager 2. The availability of spacecraft 
measurements makes the El-ring grains excellent candidates for testing aspects of the 
charging of small grains in a space-plasma environments. This is especially relevant 
since in the near future the Cassini spacecraft will make many passes through the 
E-ring region measuring the plasma energy and composition and the dust particle 
masses and velocities. 

The potential of the E-ring grains has been estimated and used as a parameter 
in recently proposed models for evolution of the E-ring particles by Horanyi et al. 
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(1992), Morfill et al. (1993) and Hamilton and Burns (1994). Horanyi et al. (1992) 
calculated the motion of the grains launched from Enceladus (a moon of Saturn which 
is the presumed source of the E-ring material) in the presence of a gravity field, solar 
radiation pressure, and electromagnetic forces. They showed that micron-sized grains, 
which are launched over a period of time and which obtain potentials between -5.4 and 
-5.8 V, would give a grain size and spatial distribution with many of the characteristics 
of the observed E-ring. Morfill et al. (1993) calculated the effect of grain potential on 
the sputtering rate of the E-ring grains. They suggested that the surrounding plasma 
is produced and maintained by "self-sputtering” of the E-ring. That is, the sputtered 
atoms and molecules are ionized and "picked up” by a planetary magnetic field and 
accelerated to corotation energies. These ions then bombard the dust resulting in a 
self-sustained process. Since the secondary electron emission coefficients, which play 
a crucial role in the grain charging, were not known for water ice, each of these groups 
of researchers made a parameter study of the effect of interest using a "best guess” 
for the secondary electron yield and the grain potential. In describing grain erosion 
Morfill et al. (1993) used a maximum yield, 6 m 1 , at energy E m as 1000 eV for the 
secondary electron emission parameters for energetic electrons incident on water-ice. 
Applying the grain-charging model described by Draine and Salpeter (1979) this gave 
E-ring grain potentials varying from -40 V at a distance 5 R, from Saturn to 1 V 
at 9 R s from Saturn, where R 3 is the radius of Saturn. Horanyi et al. (1992) used 
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the procedure described by Whipple (1981) and three values of S m with E m = 500 
eV. They estimated grain potential at the orbit of Enceladus (~4 R a ) between —8 
and —4 V. In the orbit calculations they favored 6 m = 1.5 which gives the grain 
potentials varying from -6 V (at 4 R s ) to +3 V (at 8 R 3 ). Hamilton and Burns (1994) 
calculated the motion of charged grains in the presence of gravitational forces, Lorentz 
forces, and solar radiation pressure using -5 V as the grain potential near Enceladus 
orbit. They found for this particular potential that the E-ring could sustain itself; 
i.e., charged grains comprising the E-ring strike Enceladus at high velocity ejecting 
new material. 

Using laboratory data for secondary electron yields ( Matskevich and Mikhailova 3 
1960) we recalculate values of the E-ring grain potential. We will show that extrapo- 
lation procedures for obtaining a grain potential used previously can not be done for 
the low energy plasmas in the Saturnian and other planetary magnetospheres, hav- 
ing electron temperatures usually lower than ~ 30 eV. For low energy plasmas, not 
only ’’true” secondary emitted electrons but also ’’reflected” electrons (a distinction 
based on their respective energies) constitute emitted currents from the grain and 
significantly influence its equilibrium electrical potential. In addition, the threshold 
energy for secondary electron emission, determined by the binding energy of the va- 
lence electrons, should be taken into account whenever a significant portion of the 
electron plasma is below this threshold energy: i. e., does not produce secondary 
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electrons. In section 2 we review the charging mechanism for a single-grain model 
surrounded by a plasma with Maxwellian energy distribution commonly used to cal- 
culate grain potentials. We emphasize problems that occur when models derived for 
high energy plasmas (hundreds of eV and higher) are extrapolated to the low energy 
regime. In section 3 we discuss secondary electron emission and derive the relation 
for the secondary electron current appropriate for small incident energies. In section 
4 we discuss the importance of the reflection coefficient for elastically and inelasti- 
cally reflected primary electrons, especially for low energy (cold) electrons below 20 
eV. We then obtain the charging contribution due to the reflected current and apply 
the results to Saturn’s E-ring examining the relative importance of the contibutions 
to the charging current to find equilibrium grain potentials. The plasma parameters 
based on the Voyager measurements are given in section 5. Results with applications 
to the other existing E-ring evolution models are given in section 6 and 7. 

2. CHARGING MECHANISM 

The potential of a grain depends on the energy distribution of the surrounding 
plasma as well as on the properties of the grain itself and it is determined by a balance 
between various charging currents. Assuming that the main charging mechanism 
comes from the fluxes of the incoming electrons and ions, the equilibrium potential 
of the grain is obtained from the current balance equation 
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*A "t” */sec ”1” Jrej ”1" J hv — 0 


( 1 ) 


where J e and J{ are electron and ion fluxes incident on a grain, «/,*<. is the escaping flux 
of secondary electrons, J Te { the outgoing flux of reflected primary electrons and Jh v the 
photo electron current. Distinction betwen secondary emitted electrons and reflected 
primary electrons used for calculating the grain potential has been made based on 
electron energy measurements as will be discussed. The secondary electron energy 
distribution shows a peak around E scc = 3 eV almost independent of the incident 
electron energy, while a peak which occurs around the incident electron energy is 
attributed to elastically reflected electrons. 

The secondary electron current induced by the electron impact usually determines 
the charging of a grain. Assuming spherical grains immersed in a plasma with Debye 
screening length much larger than the grain’s radius and that the fluxes of incoming 
and escaping particles are orbital-motion limited ( Laframboise and Parker , 1973), the 
incoming fluxes of electrons (e) and ions (i) have the form 


'e, = / 

Jxr 


m&x(0,±cv>) 


11 ±'- 1 )] ii iE 


(2) 


where (+ ) corresponds to electrons (ions), E is the incident electron (ion) kinetic 


6 



energy, <p is the grain’s surface potential and the factor [1 ± accounts for 

the change in the geometrical cross-section due to Coulomb attraction and repulsion 
(Spitzer, 1941). Although the electron energy distribution function in the planetary 
environments is closer to a kappa distribution ( Rosenberg and Mtndis, 1992), it can 
be reasonably fitted using two Maxwellians ( Sittler et al, 1983). Consequently, for 
calculating Saturn’s E-ring grain potential we evaluate Eq. (2) using two Maxwellian 
components for the incident electron current: a thermal (’’cold”) component with 
temperature T c and a suprathermal (”hot”) component with temperature T 

Using a Maxwellian flux distribution for incident electrons Eqs. (2) yields the well 
known result 


(J 0 exp(^m , I><01 

■■Uli + ft)/* ’ r W>o| 


( 3 ) 


with 


7 _ 1 _ 

Jq n e v e , 

4 


v, = 


! 8kT e 

7 r m e 


( 4 ) 


where n e> ,- is the electron or ion density and i > Sit - is the average thermal velocity. 

For bodies with < 0 the secondary emitted flux is obtained based on Eq. (2) as: 


Js'c {v<0) = £ + 2% *(E + e<P)dE 


( 5 ) 
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where 8(E) is the secondary electron yield. On the other hand, when tp > 0 not all 
secondary electrons will escape. Assuming that the velocity distribution of secondary 
electrons can be approximated by a Maxwellian with peak energies E iec = kT^ (2 
- 5 eV), regardless of the form of the incident velocity distribution, Prokopenko and 
Laframboise (1980) found that the escaping secondary electron flux is 

^ ( „ >0) =(1 + $-) e*p(-£-) jf <1 + f )^E + e„)dE (6) 


where the factor (1 + ^-) exp(-^^-) in front of the integral represents the fraction 
of the total electron flux emitted from the grain surface which is able to overcome the 
grain’s potential. 

While this charging scheme works for incident electron energies of the order of 
hundreds of eV, it is not appropriate when energies approach tens of eV. First, at 
low incident energies, elastic and inelastic reflection becomes the dominant process 
governing the electron loss from a grain. Also, secondary electron emission starts at 
some threshold energy (usually between 5 and 10 eV for insulators), not at zero energy 
which is often assumed in relations for the secondary electron yield 6(E) derived from 
the measurements at higher energies. Finally, in the low energy regime the escaping 
electron flux cannot be approximated by a Maxwellian independent of the incident 
electron energy because that assumption implies that for small incident energies a 
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significant portion of the secondary electrons escape with energies greater that the 
incident energy. Therefore the charging calculation needs to be modified, when the 
electron temperature is tens of eV or less. 

Secondary electron ejection can also be induced by ion or photon impact. The 
measured ion fluxes are substantially smaller than the electron fluxes and the ion 
current contribution to the charging is much smaller than the electron contribution. 
In addition, the ion induced secondary electron yield below 1 keV is small, so the 
secondary electron current induced by ion impact can generally be ignored. The 
secondary electron yield is expected to be smaller than 0.1 for ion energies in the 
order of a hundred eV ( Whipple , 1981). Therefore, the main contribution to charging 
is the direct ion current. 

For a non stationary grain Eq. (2) should be modified if the grain velocity is 
comparable to the plasma velocity ( Whipple, 1981; Havens et ai, 1987). This is the 
case for ions whose measured velocities in the inner edge of the El-ring are comparable 
to the grain’s Keplerian velocity. For the cold ion current we note that the corotating 
component of velocity is substantially larger than the thermal velocity and the average 
ion flux to a grain is roughly 


Ji « (1 — 


2 et P \ n i / v 

75) — (Vi-V kep ) 


rrii(vi - v ktp ) 2> 4 


( 7 ) 


where u,- is measured corrotating component of ion velocity and u/t ep is the grain speed. 


We will show later that the ion current contribution to the total current to a grain is 
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substantially smaller than the electron contribution. 


The photoelectron flux is also low. We use the relation given by ( Wallis and Hassan, 
1983) 


= 7 l ° u ^r ex P(- 

r AU 


max[e(/>, 0] 


) 


( 8 ) 


4 ’ r\ n ~' rv 1.3 eV 
to estimate it with the photoelectric efficiency x = 0*1 for icy grains and Saturn’s 
distance r^u =9.6. 


3. SECONDARY ELECTRON EMISSION 

Based on the energy spectra of the emitted electrons, the total electron yield, the 
mean number of ejected electrons per incident electron, is often written as a sum 

& = r + t] + 8 ( 9 ) 

of elastically (r) and inelastically ( 77 ) backscattered primaries (sometimes added to- 
gether as ’’reflected” primaries R = r + tj) and ”true” secondary electrons (6). For 
incident electrons with energies in the range where secondary electron emission dom- 
inates (for example from approximately 100 eV to few keV for ice, Fig. 2) elastically 
and inelastically backscattered primaries constitute only a small fraction of the total 
yield. In that case the total yield, <7, is often approximated by the secondary electron 
yield, S. But for small incident energies elastically and inelastically backscattered 
primaries constitute a dominant fraction of all outgoing electrons and in that case 
reflected electrons can not be neglected. 
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The secondary electron yield curve 6(E) is often described by a ’’universal” shape 
characterized by two parameters: the maximum yield 6 m and the energy at which it 
occurs, E m . Typically, the maximum electron yield 6 m is greater for insulators and 
semiconductors (1-10) than for metals (0.5-2). In the energy range where the total 
electron yield of a material is greater than one, the electron current (like an incoming 
ion current) contributes to positive charging. 

Here we use the empirical relation for the dependence of the secondary electron 
yield given by Draine and Salpeter (1979) 

S(E> = Sm <“» 

This relation approximates the secondary electron yield for normal incidence. Since 
the grain is surrounded by plasma with isotropic flux distribution an angle averaged 
yield is needed. Based on expressions given by Dionne (1973, 1975), Katz et al 
(1977) derived an angular dependence based on the range and energy loss rate for 
penetrating electrons which reads 

6(E,cosO) = 2.54 S m (£- ) L~ ex P(~ 9 . cos 0 ) (ii) 

\E m J Q cos 9 v * 

where Q is the angle of incidence of the electron and Q = 2.28 (E / E m ) 135 . This ex- 
pression gives an approximation for energies below 4 E m and for an isotropic primary 
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electron distribution it may be integrated to give an angle averaged yield ( Whipple, 


1981): 


1(E) = SMS m (4-)9~ 1+eXV( - Q) 
\E m J <3* 


( 12 ) 


3.1. Secondary electron flux for low incident energies 

The secondary electron energy distribution emitted from a material is almost in- 
dependent of the incident energy E for energies above a hundred eV and is often 
approximated by a Maxwellian with a temperature of about 3 eV ( Prokopenko and 
Laframboise, 1980) which is in reasonable agreement with measurements ( Murashov 
et a/., 1991). This is done for mathematical convinience and does not imply a thermal 
origin of secondary electron emmision. For incident energies of the order of tens of eV 
or less, we assume that it still can be approximated as a Maxwellian with a peak at 
£ S ec = £T seo but we require a cut-off at the incident electron energy. Therefore the 
velocity distribution of the escaping secondary electrons at the surface of a charged 
grain, with the required cut-off at the incident electron velocity u, is 

f(v,v,) sec = c(v) exp(-^^~) H + «p) - t;,) (13) 
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where v , is the velocity of the emitted secondary electrons at the surface of a grain, 
H is Heaviside step function and c(v) a normalization constant determined at ip = 0. 


For grains at negative potentials all secondary electrons can escape, but for positive 
grain potentials only those electrons for which \ m v] — eip> 0 can escape. In order to 
calculate J 3ec for a positively charged grain, we find the ratio A (E,<p) between escape 
fluxes from a grain at a potential ip and a noncharged grain (y> = 0) 

ME ) - -W/H’/a f ( v * v »)~ c v.-n d 3 v, 

' V So f(v,v s ) sec v s • n <Pv 3 


(» + rifc)<«p(-;fe)-(i + fg)exp(-f£*) 
1 - 


where v a • n is the secondary electron velocity component in the outward normal 
direction. This ratio represents a fraction of secondary electrons which are able to 
overcome the grain potential, i. e. escape from a grain. The total secondary electron 
flux in for ip > 0 in Eq. (6) becomes 

Jsec,(v>o ) = J Q (1 + -jjr) S(E + eip)dE (14) 

When the incident energy is substantially larger than the grain potential, the factor 
reduces to (1 + k ) exp(— which can be pulled in front of the integral 
so that Eq. (14) becomes equivalent to the result in Eq. (6). 
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3.2. Secondary electron emission for H 2 0 - ice 

Unfortunately, very few measurements have been made of the secondary electron 
yield for water-ice, or for other molecular ices of interest in the outer solar system, 
especially for low incident energies. Therefore, modellers have used different param- 
eters, 6 m and E m , as well as different energy dependencies, S(E) 1 inferred from the 
yields measured for other materials. 

Recently, Suszcynsky et ai (1992) measured the secondary electron yield of an 
H 2 O ice film at normal incidence using a scanning electron microscope. Incident 
electron energies were between 2-30 keV, far above the maximum-yield energy E m . 
Suszcynsky et ai (1992) used the Sternglass universal curve to extrapolate measured 
values in the low energy range, predicting the secondary electron emission parameters 
for water-ice of 8 m =6.8 and E m =142 eV. These authors were apparently unaware of 
an earlier measurement of the secondary emission yield by Matskevich and Mikhailova 
(1960) for electron energies from 100 - 2500 eV at normal incidence, using a single 
pulse method. 

Both measurements are plotted in Fig. 1 together with the numerical fit to the data 
using functional dependences given by Eq. (11) with 6 = 0 (curve a) and (10) (curve 
b), for the secondary electron coefficients 6 m =2.35 and E m =340 eV. The widely 
used Sternglass furmula [Meyer- Vernet, 1982) is also plotted in the same figure (c) 
for comparison. Above 1 keV the different fitting formulas give significantly different 
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yields with Eq. (10) giving the best approximation to the measurements. For the low 
energy portion of the secondary electron yield curve, which is most important in the 
E-ring, the fitting formulas do not differ greatly, but the data is very different from 
the ’’best guesses” of 6(E) used by Morfill et al. (1993) and Horanyi et al. (1992) 
to calculate the grain potential (Fig. 1 ). We use Eq. ( 11 ) and, as a check Eq. ( 10 ), 
both used in recent models ( Morfill et al. 1993, Horanyi et al. 1992) to see how much 
the calculated grain potential is affected by the different functional relationships for 
6(E). 

As the ’’temperature” of the secondary emitted electrons, we use the peak energy 
~3 eV of the secondary-electrons emitted from quartz (which has secondary electron 
parameters 6 m and E m similar to mica, glass and water-ice), measured by Murashov 
et al. (1991). 

3.3. Threshold energy for secondary electron production 

When dealing with low energy electrons inside the solid it is important to notice 
that only those absorbing sufficient energy from an incoming electron, ion or photon 
can leave the surface and contribute to the secondary current. The surface barrier 
for insulators is determined by the electron affinity (E J 4 ) which is the energy differ- 
ence between the vacuum level and the bottom of the conduction band. Only those 
electrons for which the component of kinetic energy perpendicular to the surface is 
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greater than Ea will escape from the material. In addition, electrons from the va- 
lence band need to absorb at least the band gap energy, which is 7.8 eV for cubic ice. 
Since the electron affinity is 0.9 eV, the valence band edge lies about 8.7 eV below 
the vacuum level for cubic ice ( Baron et al, 1977) or about 9 eV for amorphous ice 
( Williams et al, 1974). The primary electron also gains ~0.9 eV (electron affinity) 
when entering material, so we put the threshold energy for secondary electron pro- 
duction at E t h = 8 eV. In Figure 2 we plot the functions in both Eqs. (10) and (11), 
starting at E t h = 8 eV which corresponds to replacing E with E — E t h and E m with 
E m — E t h in Eqs. (10) and (11). 

A lower threshold may result for the excitation of electrons in trap states in the 
band gap with binding energies ~ 1 — 2 eV ( Haas et al , 1983). If the grain is negatively 
charged, then, of course, the excess electrons must reside in traps, which are likely 
to be near the surface. Although these electrons can be more easily removed from 
the grain, they constitute an extremely small fraction of the electrons involved in 
determining the current balance. That is, for a 1 micron radius grain at -10 V there 
are only 5 x 10 5 excess electrons per surface molecule. Therefore, we will use the 
threshold value given above for all secondary electrons. 
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3.4. Secondary electron emission for isotropic incidence 


The measurements described above were for normal incidence. Morfill et al. (1993) 
used a multiplying factor of 2 for secondary electron yield to account for both the as- 
sumed spherical shape of a dust particle and isotropic incidence ( Draine and Salpeter , 
1979). Measurements show [e. g., Salehi and Flinn, (1981)] that the enhancement of 
the secondary electron yield due to isotropic incidence has a different energy depen- 
dence and cannot be accounted for by a simple multiplying factor. It was observed 
that the secondary electron emission increases with incident angle and that the value 
of E m shifts toward higher energies, approximately proportional to (cos#) -1 . The 
enhancements of the secondary electron yields due to the small particle effect ( Chow 
et al., 1993), which can be significant for a grain size of order 0.1 pm and smaller, 
can be ignored in the case of the E-ring grains. 

First we use Eq. (11) as a fit for the normal incidence yield and find the angle 
averaged yield from Eq. (12). In Fig. 3 the dash-dot lines represent normal incidence 
yield from Eq. (11) (lower line) and the corrected yield for isotropic incidence given 
by Eq. (12) (upper line). As a check we also find an angle averaged yield using Eq. 
(10) as the normal incidence yield. For this we use angular dependent measurements 
( Salehi and Flinn, 1981) for V 2 Os — PiOs — Cs 2 0 glass which has secondary emis- 
sion parameters close to ice. We numerically integrated those measurements over all 
incident angles and find an angle averaged yield which we can scale to ice, assuming 
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that parameters E m and S m for ice and glass in the case of the isotropic incidence are 
shifted by the same factor from the normal incidence parameters. We plot in Figure 
3 as solid lines the normal incidence yield given by Eq. (10) (lower line) and the 
calculated angle averaged yield based on Eq. (10) (upper line). 

4. REFLECTION COEFFICIENT 

A typical ejected electron flux distribution, measured by Harrower (1956) for tung- 
sten target at the incident electron energies of 10 and 20 eV, is shown in Fig. 4 (solid 
lines). The first peak at ~ 3 eV corresponds to secondary electrons while the second 
peak at the incident energy corresponds to the reflected electrons. The maximum is 
elastic reflection while inelastic reflected electrons, which lose some of their energy in 
the interaction with the target, correspond to energies below the elastic maximum. To 
model such flux distributions at other incident electron energies we use two curves: 
one for the secondary electrons and the other for reflected electrons and calculate 
these two current contribution (J sec and J rc j ) independently. The distinction be- 
tween the inelastically scattered electrons and the secondary electrons is somewhat 
arbitrary, but both contributions, J sec and J rc /, enter into the total charging current 
as a sum. The outgoing energy only determines whether the electron can escape from 
a grain’s potential well, i. e. how it contributes to the charging. 

For negative grain potentials < 0) all emitted electrons can escape regardless 
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of their energy and the reflected flux in that case can be calculated using Eq. (5), 
where S is replaced by tj + r — R, and where we add both elastically and inelastically 
reflected electrons in the reflection coefficient. 

For positive grain potentials only electrons with sufficient energies escape, i.e., all 
elastically reflected electrons and those inelastically (as same as secondaries) emitted 
with energies v? > tp. To see how the energy distribution of reflected electrons 
influences the current to a grain, and consequently the grain potential, we present a 
simple model for reflected current from a charged grain. For the sake of simplicity, 
we approximate all reflected electrons from an uncharged grain with a half-Gaussian 
velocity distribution up to the elastic peak at incident electron velocity v 

frc}(v , U r ) = C(v) exp J^ ) H{v - V T ) (15) 

where v T is velocity of the reflected electrons at the surface, T rc j measures the spread 
of the distribution around the elastic peak, while H{y — u r ) introduces the cut-off 
at the incident electron energy and c (v) is normalization constant. We also assume 
that the electrons are reflected from the small surface element isotropically ( Whipple 
and Parker 1969, Whipple 1981), which means that reflected 'fluxes measured from 
the surface normal follows the experimentally known "cosine law". We use the same 
procedure as that to obtain Eq. (14) for the secondary electron emission flux to 
estimate reflected electron contribution J re j in Eq. (1). The total reflected flux for 
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where R is a measured reflection coefiicient. The dashed lines in Fig. 4 represent our 
approximation for the flux distribution which includes secondary electron flux (with 
peak at kT^ = 3 eV) and reflected flux with our choice for kT re j = 1 eV. Due to the 
finite resolution of the energy analyzer the measured energy spreads around the inci- 
dent energies (10 and 20 eV) are overestimated ( Harrower , 1956), and, consequently, 
our choice for the parameter kT re f represents the upper limit for the actual spread 
around the elastic peak. However, the change in the equilibrium grain potential in- 
troduced by this choice does not exceed 0.1 V, which is less than other estimated 
uncertainties. It can be concluded that when the elastic peak is sharp (which is the 
case for water ice, Michaud and Sanche , 1987b), the shape of the reflected electron 
energy distribution does not significantly influence the charge balance. 

4.1. Electron reflection from H 2 O ice 

From measurements of the reflection coefficient for various materials at low incident 
electron energies ( Bronshtein and Novitskii , 1978; Nemchenok et al, 1976; Khan et 
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a/, 1963; Fridrikov and ShuFman } 1959) it is known that the reflection coefficient 
approaches zero as E — * 0, reaches a maximum below "20 eV, and decreases slowly at 
higher energies. Typically, for metals the maximum value reaches 0.1 - 0.4, whereas 
for dielectrics it attains 0.5 - 0.8 ( Dobretsov and Gotnoyunova, 1971). 

Measurements of the elastic electron reflectivity of amorphous films of water-ice 
have been done by Michaud and Sanche (1987a, 1987b) for various thicknesses (1-40 
monolayers) and for incident electron energies from 1 to 20 eV at 14° incidence. Using 
the same experimental arrangement as Michaud and Sanche (1987a, 1987b), Bader et 
aL (1988) measured transmitted current on a 50-layer H 2 O film for incident electron 
energies from 0.1 to 4 eV. Matskevich and Mikhailova (1960) measured reflection 
coefficients for ice in the energy range from 100 eV to 2.5 keV. In Figure 2 we show 
both measurements together with a least squares fit to the measurements (dashed 
line) using, for convenience, a form similar to that in Eq. (10), which we use as R(E) 
in calculations. This fit should be considered only as a rough approximation to the 
actual functional dependence R(E) in the range of interest, below 1 keV. 

5. PLASMA PARAMETERS 

For equatorial plasma parameters around Saturn’s E-ring we use the model given 
by Richardson and Sittler (1990), which is based on Voyager measurements. The elec- 
tron energy distribution is fitted by two Maxwellian distributions: thermal (’’cold”) 
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component with temperature T c and suprathermal (”hot”) component with temper- 
ature Th- The high energy tail is not well fit using two Maxwellians, but it does not 
play an important role in our calculation, since the electron fluxes and the secondary 
electron yields are low in this region. For the proton and oxygen ion densities, we use 
values extrapolated to the equatorial plane using the Richardson-Sittler model (1990). 
A Summary of the plasma parameters used are given in Table 1, where T ec , T e ^, n ec , 
and n e fc represent temperatures and densities of the cold and hot electrons respec- 
tively. There are, of course, considerable uncertainties in this model, which involved 
extrapolation of the measured plasma densities to the equatorial plane, and indeed 
the plasma parameters probably vary in time somewhat. Both the ion and electron 
temperatures recorded during the Voyager 1 and 2 traverses vary substantially along 
the same dipole L-shell. Consequently, the distribution of temperatures and densi- 
ties can be considered a parameter study for grain charging since the distance from 
Saturn affects the result only via changes in electron temperature and density. Since 
Voyager 1 crossed the equatorial plane at L ~ 6 Saturnian radii, equatorial plasma 
parameters can be determined there and are also given in Table 1. 

6. RESULTS: POTENTIAL OF THE E-RING GRAINS 

Using Eqs. (8), (3) and (7) for the photon, electron and ion currents, Eqs. (5), (14) 
for the secondary electron currents and Eqs. (5), (16) for the reflected current, a grain 
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potential can be found for which the current balance equation (1) is satisfied. As an 
example, in Figure 5 we give all currents to a micron-sized ice grain at radial distance 
L = 6 R s , incident electron current to a grain, secondary electron current, reflected 
current, sum of the photon and ion currents for H + and 0 + ions and, finally, total 
current to a grain. Current balance is obtained in this case for the grain potential 
-2 V. As one can see the combined photoelectron curent and the ion current (due to 
the small ion velocities) to a grain is substantially smaller than other contributing 
currents; consequently, the equilibrium potential primarily results from balancing 
different electron currents. 

In Figure 6 we show the calculated potential for a spherical grain as a function 
of radial distance from Saturn for both relationships used for the secondary electron 
emission (a, b in Fig. 1 and Fig 2). The resulting potential using Eq. (11) for normal 
incidence yield (and, consenquently, Eq. (12) for the isotropic yield) gives an average 
''-'IV more negative potential than that based on Eq. (10), since the secondary 
emission yield in that case rises more rapidly up to E m (Fig. 3). The potential in 
Fig. 6 is seen to rise from 4-8 R a and stays almost constant between 8 and 10 R a . The 
potential increases in going from Enceladus (4 R s ) to 8 R s due to increasing electron 
temperatures (both hot and cold) and increasing densities of the hot electrons. The 
almost constant potential near Rhea (8.7 R s ) and further outside arises from the 
fact that, when the positive potential becomes greater than the secondary electron 
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temperature, JcT scc = 3 eV, the number of secondaries (which constitutes the major 
part of the positive current to the grain) sharply decreases. These results are also 
compared to the best guess of Morfill et al. (1993) and three models from Horanyi 
et al. (1992). Agreement is seen to be best with model H2 in which S m = 2. In their 
calculation of E-ring grain orbits Horanyi et al (1992) used S m =1.5 and E m =500 
eV, whereas Morfill et al. (1993) used S m =1 and E m =1000 eV to describe E-ring 
grain erosion. The best fit of the experimental yield is obtained with 6 m =2.35 and 
E m =340 eV as shown in Fig. 1. 

Previous authors have not separated the reflected current in calculating the charg- 
ing or in calculating the equilibrium potential, which from Fig. 5 is clearly seen to 
be important at almost all potentials shown. Not including reflection, the E-ring po- 
tential would appear 3—5 V more negative than that calculated throughout the full 
range of distances. Since there are no measurements of the secondary electron yield 
or the reflection coefficient for ice for a few eV to tens eV, the estimated reflection 
coefficient in that range (discussed above and given in Fig. 2) is a source of uncer- 
tainty for the calculated grain potential. The threshold energy for secondary electron 
emission also impacts calculated potential when the plasma temperature is low. For 
instance, not including the threshold energy for the secondary electron production 
(we use 8 eV) would give potentials 0.5 — 1 V more positive than those calculated. 
Since most of the inelastically reflected electrons from ice at low incident energies 
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suffer small energy losses, all reflected electrons can be considered as having reflected 
elastically. The choice, kT re j = 1 eV (as opposed to kT re j = 0 eV, if all electrons were 
reflected elastically), introduces a change in the grain potential of the order of 0.1 V. 
Further uncertainties may result from the effects of irregular shapes and roughness of 
ice grains, impurities in the ice and radiation induced defects which introduce electron 
states in the band gap, thereby affecting the yields and threshold energies. Using the 
Voyager 1 ring plane crossing data only (Table 1) the grain potential at 6 R s would 
be +3 V as compared to -2 V for value b in Fig. 6. Therefore, the plasma parameters 
are the largest uncertainty in this calculation. 

7. CONCLUSION 

In this paper we calculated the charging of an ice grain in Saturn’s magnetosphere 
where the plasma electron temperatures are low. These calculations were based on ex- 
trapolations of the complete set of available secondary electron and reflected electron 
data, taking into account the physics of the secondary electron yields and electron 
reflection coefficients. We assumed that the grains are one micron solid water-ice 
spheres, as suggested by most of the measurements ( Showalter et. a/., 1991.). How- 
ever, the calculated potentials are not sensitive to the grain size as long as the grains 
are not much smaller than ~ 0.1 y.m and the secondary electron yield parameters for 
other materials possibly present in the E-ring probably do not differ significantly from 
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that of water ice. We show that for low incident electrons energies the calculation of 
the grain potential requires consideration of electron reflection and of the secondary 
electron production threshold. In fact, for the potentials calculated here for Saturn’s 
inner magnetosphere, varying from about -5.5 V at 4 R s to 5 V at 10 R a , reflection is 
always important, as it is in the case whenever the electron temperature is < 20 eV 
or less. Therefore, the charging calculations developed here can be used when the 
plasma parameters in Saturn’s inner magnetosphere are more firmly established and 
in calculations for other planetary plasmas environments in which a cold electron 
component is dominant. 

Our calculation of the E-ring grain potential vs. distance from Saturn using the 
plasma data of Richardson and Sittler (1990), which is an average of Voyager data with 
ion densities extrapolated to the equatorial plane, is seen to be in rough agreement 
with curve H2 (Fig. 6) in Horanyi et al . (1992) who described the spatial distribution 
of the E-ring. Therefore, our results appear to support this aspect of the hypothesis 
of the formation and evolution of the E-ring grains. On the other hand, our results 
differ significantly from the best guess of Morfill et al. (1993) for describing a plasma 
source by low energy ion sputtering of E-ring grains. However, sputtering of the E- 
ring grains by keV ions does contribute significantly to the plasma formation near 
Enceladus ( Johnson et a/., 1993). 

Horanyi et al. (1992) and Hamilton and Burns (1994) concluded that the spatial 
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distribution and grain size characteristics of the E-ring could be understood by grains 
being launched from Enceladus, becoming charged in the ambient plasma, and or- 
biting under gravitational and electro-magnetic forces. Horanyi et al. (1992) showed 
that if micron-sized grains launched over a period of time from Enceladus (3.95 R ,) 
would obtain potentials ~ —5.4 to —5.8 V they would disperse over time producing 
an optical depth profile with a thickness like that of the actual E-ring. In spite of 
the fact that their "best guess” estimate of the secondary electron coefficients, made 
in order to obtain the observed E-ring characteristic, underestimates the secondary 
electron yield below 1 keV, agreement with our potentials occurs because they also 
neglected the electron reflection coefficient. Hamilton and Burns (1994) found that 
micron-sized grains charged to -5 V would obtain the required orbital eccentricity to 
account for the spatial extent of the E-ring. That is, they would cross the orbit of 
Tethys during their life in Saturn’s inner magnetosphere. Such potentials are close to 
those found here at 4 R a , but the fact that the potential changes with distance from 
Enceladus must be considered. Therefore calculations presented here appear to con- 
firm aspects of the E-ring hypothesis, if the plasma parameters used are reasonable. 
These results can now be used for more detailed determination of the physics of the 
E-ring. 
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FIGURE CAPTIONS 


Figure 1. Secondary electron yield for ice vs. incident electron en- 
ergy (keV): measured by Matskevich and Mikhailova (x’s) and 
Suszcynsky et al. (diamonds); functional dependences given by 
Eq.(ll) with 6 = 0 (a), Eq.(10) (b) and by the ”Sternglass law” 
8(E) = 7.4<5 m ^ r exp(-2 v /^) (c) with 8 m =2.35 and E m =340 
eV; yields used as the ’’best guesses” by Morfill et al. (M) and 
Horanyi et al. (H). 

Figure 2. Measured secondary electron yield given by Matskevich 
and Mikhailova (x’s) and fitted yields from Eq. (10) (b) and Eq. 
(11) with 6 = 0 (a), both starting at E t h = 8 eV as the threshold 
energy. The reflected yield measured by Bader et al. (squares) and 
Matskevich and Mikhailova (-f-’s) is shown together with our fit to 
the actual functional dependence R(E) = (1 (dashed 

line). 

Figure 3. Secondary electron yield for ice for normal incidence from 
Eq. (11) (lower dash-dot line) and the corrected yield for isotropic 
incidence given by Eq. (12) (upper dash-dot line). Solid lines 
represent the normal-incidence yield given by Eq. (10) (lower line) 
and the calculated angle averaged yield based on measurements on 
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glass (upper line). 

Figure 4. Total electron fluxes from tungsten for 10 and 20 eV pri- 
mary electrons at normal incidence measured by Harrowtr (solid 
lines) and modeled flux distribution (dashed lines) using velocity 
distributions in Eq. (13) for secondary electrons and that in Eq. 
(15) for reflected electrons. 

Figure 5. All currents to a micron-sized ice grain at radial distance 
L = 6 Saturnian radii: electron current incident to a grain (I e i), 
secondary electron current (I , wc ), reflected electron current (/ r «/), 
sum of the photoelectron and ion currents for H + and 0 + ions 
(f»on+pAot) 5 and total current to a grain (I tot)- The photoelectron 
and the ion currents are both small, of the order of 10 -19 A. Cur- 
rent balance is obtained for the grain potential —2 V. 

Figure 6. Calculated grain potential as a function of the radial dis- 
tance in Saturnian radii. Solid lines show the potential determined 
in this paper using E m =340 eV and 8 m =2.35. Two different 
extrapolations for the secondary electron yield lead to slightly dif- 
ferent potentials: line (a) based on Eq. (11) and line (b) based on 
(10). The dash-dot line represents the potential favored by Mor- 
fill et. al. (1993) with estimated secondary electron parameters 
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E m =1000 eV and S m =1. Dashed lines show potentials given by 
Horanyi et. al. ( 1 992) using E m =500 eV and three different val- 
ues for the maximum yield: 6 m =0 (HO), S m =1 (HI), and S m =2 
(H2): their ’’best guess”, based on the observed E-ring character- 
istics, is 6 m =1.5. 
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Table 1. Equatorial Plasma Parameters 1,2 


L 

T 

A cc 

n ec 

T e fc 

H eh 

T 0+ 

n 0 + 

T„ + 


4 

2.7 

90 

100 

0.2 

40 

70 

12 

20 

5 

3.5 

45 

120 

0.4 

80 

40 

14 

3.7 

6 

5.0 

27 

150 

0.4 

100 

25 

16 

2 

7 

6.8 

15 

170 

0.4 

120 

15 

18 

1 

8 

11 

4.5 

200 

0.4 

170 

4 

20 

0.8 

9 

13.5 

2.7 

250 

0.3 

220 

2.3 

22 

0.65 

10 

17 

1.9 

300 

0.2 

260 

1.6 

24 

0.5 


Table 1. 1 Plasma parameters from Richardson and Sittler (1990): 
temperatures T (eP) and densities n (cm -3 ) for electrons (”cold” 
and "hot”) and ions ( 0 + and H + ) in equatorial plane versus 
distance in Saturnian radii L 

To indicate uncertainties in plasma parameters we used 
Voyager 1 measurements ( Sittler et al ,1983) averaged over a two 
huor period covering the ring plane crossing between 5 < L < 
6.7. The proper electron density is determined from ion density 
requiring the charge balance. These parameters are: T ec = 12.3 ± 
2.9, n ec = 24.3 ± 6.2, T efl = 99 ± 50, n eh = 1.0 ± 0.4 
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Particle Emission Induced by Ionization Tracks in Water Ice 
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abstract 


We study electron emission and sputtering during irradiation of amorphous water ice at 60 K by 

, Li , Be , B , C , N*, O , F , and Ne* ions in the energy range from 10 to 100 keV. 
We find that for constant velocity (5 keV/amu) ions the dependence of the sputtering yields with 
projectile atomic number 2 is proportional to the square of the Z-dependenoe of the electronic 


stopping power dE/dx, using dE/dx values predicted by Yarlagadd* Robinson, and Brand,. The 
electron yields increase sublineariy with dE/dx, indicating the strong influence ofun-neutralized 
holes in the ionization track produced by the projectiles. 


•Pennanen, address: OFD, Rise National Ubomoty, Denmark 
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1. INTRODUCTION 

Fast ions moving through condensed matter produce a track of excitations and ionizations 
that can have lasting effects depending on the properties of the medium [1, 2, 3, 4]. A core of 
unbalanced holes results [5] from the fact that electrons receive almost all of the energy transfer 1 

which ejects them from the path of the ion. The central problem in describing the track is to 
understand the evolution of the electrons, holes, and excitations with time after the primary 
excitation events. In a metal the ionizations are screened very quickly by the free electrons (within 
10 ,6 -10 15 s), whereas the ion core in an insulator may survive for much longer times (>10' n s). If 
the density of holes and liberated electrons is high, the behavior of the charged particles in the 
track becomes extremely complex, until neutralization eventually occurs. 

The persistence of an unbalance of excess holes in the core of a track in insulators can 
damage the material through the action of the very large Coulomb repulsion between holes. This 
type of radiation damage is important in radiation biology, where it may cause the destruction of 
living cells [2,3], In certain materials the damage tracks can be made visible and counted, and the 
information used in radiation dosimetry and in the dating of rocks [1], 

The emission of particles (electrons, ions, neutral species) from the track of a projectile 
entering the insulator [6, 7, 8, 9] depends on ionization events that take place not only 
immediately below the surface, but further into the material as well. Because ejected electrons 
cannot return to the track to recombine with the holes, the effects of unbalanced Coulomb forces 
persist until the holes diffuse away (usually slowly) eventually reaching a conducting electrode 
where they neutralize. 
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The important process of ejection of electrons from insulators is not well understood due 
to uncertainties in the intervening physical mechanisms. Electrons may be excited into the 
conduction band by transitions from the valence band or core levels either produced directly by 
the projectile or through cascade multiplication of the excited electrons [10, 1 1] As is common 
to other ionization processes, the energy distribution of the excited electrons falls off very rapidly 
with electron energy [12], Most of the excited electrons will have insufficient energy to excite 
other electrons across the band gap so they will lose energy slowly through excitation of atomic 
vibrations as they move away from the point of ionization. For dilute tracks, this motion will be 
mainly affected by the residual parent hole and by strong elastic scattering with atomic cores 
which tends to randomize the electron motion. Those electrons with a sufficiently large kinetic 
energy component perpendicular to the planar surface barrier (of magnitude equal to the electron 
affinity of the surface) will be ejected from the solid. The time between excitation and emission is 
expected to be small, of the order of a few fs, the value estimated for metals [13]. In this limit of 
dilute tracks, the number of emitted electrons is proportional to the number of electrons excited 
per unit path length, which is in turn proportional to the electronic stopping power of the 
projectile (dE/dx) at high velocities [10,1 1], 

The termination of the track on the surface implies that the events which lead to radiation 
damage in the bulk of the solid can now eject atomic and ionic species, though most are neutral 
[14]. This form of sputtering can be much more important than elastic sputtering by direct 
momentum transfer by the projectile to lattice atoms, particularly for condensed gas solids [7]. 
Unlike elastic sputtering, where ejection occurs typically within a ps after the projectile hits the 
surface [15], electronic sputtering can occur over a much longer time scale (> ps). The energy 
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required for atomic motion does not necessarily come from Coulomb repulsion between holes 
(this process, in fact, has not yet been conclusively demonstrated) but from the decay of excitons 
into repulsive neutral states or, in molecular solids such as ices of water or carbon monoxide, 
from more complex processes involving radicals [16, 17, 18, 19J. 

The case of water ice is particularly interesting because of multiple applications. Among 
them are radiation biology, since amorphous ice is a good analog to liquid water, and 
astrophysics, where ices on cosmic grains and outer solar system objects are bombarded by 
ionizing radiation [19]. Brown et al. [17,20] have shown that ice can be sputtered readily by fast 
protons and found that the sputtering yield is proportional to the square of the electronic stopping 
power at MeV energies but deviates from this behavior for protons of tens of keV with similar 
values of electronic stopping power [20], 

The energy release processes that lead to particle ejection during electronic sputtering of 
ice are not known in detail. Presumably, the ionizations and some of the higher-energy excitations 
initiate a sequence of processes from which one or several radicals are generated. When radicals 
from different ionization events react, some of the processes are sufficiently exothermic that intact 
molecules or fragments can be ejected [7, 19, 21, 22,]. 

In this work we present studies of the dependence of electron yields and sputtering on 
electronic stopping power for H\ D\ He\ Li\ Be + , B\ C\ 1ST, 0 \ F + , and Ne + ions in the energy 
range from 10 to 100 keV. We analyze the data in terms of the density of the ionization track. 

2. EXPERIMENTAL 

Experiments were performed in an ultrahigh vacuum chamber (Fig. 1) connected to a 120 
keV ion accelerator. The base pressure in the chamber was ~10* 10 Torr , rising to 2 - 5 x 10' 9 
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Torr during the measurements. The rotatable target assembly is cooled by a closed-cycle 

refrigerator which can reach 20 K. A gold-coated quartz-crystal microbalance positioned at the 

target is used to measure the sputtering yields and film thickness. A second crystal mounted 

behind the target (not shown in fig. 1) is used in a heterodyne method as a reference to 

compensate the temperature dependence of the crystal frequency [23], In this system, the areal * 

mass sensitivity corresponds to about 0. 1 monolayer of ice. Surrounding the target is an 

aluminum cylinder which acts as an anode for electron collection (when biased at +300 V) as well 

as a thermal shield against the heat from the chamber walls. Between the target and the anode is a 

92 /o transparent nickel screen which is biased at -90 V from the anode to suppress secondary 

electron emission from the anode. 

Amorphous ice films of 1 200- 1 500 A were grown by flowing vapor of outgassed high 
purity water through a capillary array doser onto the target crystal held at 60 K, then bombarded 
by 10-100 keV ions at normal incidence. The ranges of these ions in ice are all greater than the 
film thickness. Sputtering yields are determined by measuring the ion beam dose (typically -10 14 
ions/cm 2 ) and the frequency change of the crystal. The electron yields (possibly including a very 
small fraction of negative ions) are determined by measuring the ion beam current to the target 
assembly and the electron current to the anode. 

3. RESULTS AND DISCUSSION 
3 A. Sputtering 

The sputtering yield Y s from water ice induced by 10-90 keV protons is shown in fig. 2. 

Sputtering induced by ion bombardment depends primarily on the binding energy of surface atoms 



6 


and the stopping power of the incident ion [7], The dependence of the sputtering yield on the 
electronic stopping power is often analyzed in terms of a power law 

r s = c(dE/dx)" ( 1 ) 

where c is a constant for a given material and n is determined empirically. The value of n is 
related to different physical processes: n= 1 implies that sputtering results primarily from 
individual excitation or ionization events, while n= 2 implies that excitations or ionizations interact 
to give rise to sputtering. Values in between and beyond indicate that varying degrees of 
interaction processes are responsible for sputtering. For our data in Fig. 2 we obtain n - 1 .3 ± 0.2, 
using the experimental dEJdx values of Bauer et al. [24] which we extrapolate below 30 keV. 

In fig. 3 data for sputtering as a function of atomic number of the projectiles are 
presented together with the corresponding stopping powers, at a constant velocity corresponding 
to 5 keV/amu (9.8* 10 7 cm/s). The square root of the sputtering yield as well as the stopping 
power have been normalized to the proton value. Since there are no measurements of stopping 
powers of ice for heavy projectiles, we use the Z-dependence of the stopping power given by 
Yarlagadda et al [25], which have been shown previously to describe accurately the Z-dependence 
of the electron yields from aluminum [26], valid in the limit of low velocities. The Z-dependence 
of stopping power (dE/dx) z for slow ions at constant velocity is given by: 

(dEjdx) z = CZ 2 |l - exp(-0.95Z‘^ )j 2 ( 2 ) 

which describes the overall trend but not the structure apparent in the experimental data of fig. 3 
as Z goes from 5 to 6. We believe this structure is due to shell effects (Z-oscillations [26]) in the 
stopping power. Unlike the case of sputtering by protons of different energy, we see that 
changing dEJdx at constant velocity produces a quadratic dependence of the sputtering yield with 
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stopping power, as found for MeV protons. These two apparently inconsistent results can be 
reconciled by a velocity dependent c factor in Eq. (1). 

3 B Electron Emission. 

It is well known from the electron emission literature on metals that the electron yield is 
proportional to the stopping power [10, 26, 27, 28, 29], This is a direct result of the fact that the 
electrons are dominantly produced by cascades rather than primary ionization alone [10, 30] 
which is proportional to the ionization cross section. Primary ionization alone is important for 
slow projectiles on insulators with a large band gap, like water ice. However, since the 
dependence of the stopping power [24] and the ionization cross section [31] with proton energies 
above 10 keV are fairly similar, and since the ionization cross section for other ions is not 
available, we will use the existing stopping power data. 

The electron yields for protons increase with energy from about 1 .7 electron/proton at 5 
keV up to about 3 electrons/proton around 70 keV in the region of the maximum of dEldx (fig. 

2). The yields increase slower with energy than expected from the variation in the stopping 
power; we will return to that point below. The electron yield was found to be independent of the 
beam current density over an order of magnitude around the current densities used in these 
measurements (~ 0.2 pA/cm 2 ) , indicating that electron emission from these thin films is a single 
incident ion effect and does not involve macroscopic charging. 

The electron yield for a number of ions incident on water ice with energy 5 keV/amu is 
shown in figs. 3 and 4 as a function of the atomic number of the projectile, Z. The yield increases 
towards a value around 3.6 electrons/ion for larger atomic numbers. The electron yield for a 
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number of ions incident on aluminum [26] is also shown in fig. 4 to compare to the case of metals. 
One notes that y for aluminum is much lower than that from water ice for light projectiles, but that 
the ratio y(ice)/y(Al) decreases to a value around 1 .4 for high Z. The lower electron emission for 
aluminum than for ice is consistent with the fact that metals generally have a low y compared with 
insulators, but the strong decrease in the yield ratio suggests that an additional effect due to high 
ionization densities is acting in the insulator. 

The Z-dependence of the electron yield from aluminum has been shown to be described 
very well by the stopping power expression from Yarlagadda et al. [25] mentioned above. 

Unlike the case of Al discussed by Alonso et al. [26], where there is a good agreement between 
the Z-dependence of the electron yields and the Z-dependence of the stopping power, for water 
ice, the y(Z) dependence lies far below the stepping power curve (fig.3). Therefore, one is led to 
the conclusion that electron emission from water is not just determined by the stopping power, 
but is strongly influenced by track processes. The decisive difference between Al and ice is that, 
unlike the case of Al, the ions in the track of ice are not screened during the electron ejection 
process. As indicated schematically in Fig. 5, the unbalanced holes in ice attract the electrons, 
preventing them from escaping, and, consequently, reduce the yield for a high hole density, i.e., a 
high stopping power. The magnitude of this effect caused by the unbalanced holes will be studied 
quantitatively in the following section on the basis of a simple model. 

3 C. Model for Track Effects in Electron Emission 


The electron yield from a metal can be well approximated by [10, 28]: 
y = BdE/dx 


( 3 ) 



B is a material dependent factor that takes into account the energy distribution of electrons in the 
solid and the transmission of electrons through the surface barrier: 
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B = \f(E)T(E)dE (4) 

u 

As mentioned in the introduction, the internal energy distribution of excited electrons, J(E), is a 
decreasing function ofE, but the precise shape is not known for ice. We model it as flE) = b/E 2 
where b is a constant and assume that this functional form is independent of the track density. 
Since low energy electrons undergo strong elastic scattering inside the solid, they will arrive at the 
surface as an isotropic flux in the half space inside the solid. The probability of transmission 
through the planar surface barrier of magnitude U is given by [10]: 



The sub-linear increase of the electron yields with dEldx observed here has been reported 
by Jacobsson and Holmen [32, 33] for ion bombarded Si02. These authors proposed that it is due 
in part to an additional barrier created by the space charge of the track. They fitted their data 
using an additional planar surface barrier to account for the space charge. We attempt here to 
calculate the average track potential acting on the escaping electrons. Furthermore, we assume 
that this potential will act similarly to a spherical barrier which is added to the planar barrier that 
the electron must overcome to escape. The spherical barrier is included by changing the limits of 
integration in Eq. (4) and by simply shifting the electron energy in the transmission function by Us. 
Then: 



t 


( 6 ) 



The form of the transmission function accounts for both barriers: for Us= 0 the transmission 
function is identical to the planar value (Eq. (5)), whereas for U, » U, it approaches the 
standard expression for a spherical barrier of height U , . Integration yields 
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( 7 ) 


which tends to y = — — — for large track potentials. 

U s dx 

Calculation of an effective potential U s is based on the assumption that the incident 
projectile produces a line of holes from the surface to the substrate separated by W/(dEJdx) where 
W is the mean energy required to create an electron-hole pair (29 eV for water [34]). The 
potential is calculated at the surface by assuming that the number of unbalanced holes is equal to y 
and that the potential is produced from the total charge placed at a depth LI 2 , where 
L = yW/( dE/dx ) . This leads to a potential Us 


v _ Iq 1 dE/dx 


Ane W 


( 8 ) 


where the dielectric constant ej e 0 —1.7 [35], appropriate for energy transfers of eVs and tenths 
of eV. Using the density of 0.82 g/cm 3 for vapor deposited water ice [36] we obtain 
Uf=0.\6\ (dE/dx) (eV/(10 l5 H2O/cm 2 )), which gives a track potential of ~9 eV for the highest 


stopping power used. 

Figure 6 shows the results of plotting y = B(dE/dx) t from eq. (7) with the data from 

figures 2 and 3. The stopping cross sections used for the proton data are from Bauer et al. [24] 
and those used for the other ions are obtained from the proton data using eq. (2). For these 
calculations we applied the potential Us obtained from eq. (8) and the surface barrier was taken to 


s 
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be U— 0.9 eV, equal to the electron affinity of water ice [37], The curve was normalized to a 
value of 1 .9 electrons/ion at 10 eV/(10 15 H2O/cm 2 ). Even though the model is based on several 
simplifying assumptions, it describes well the trend of the data over the entire stopping power 
range. We believe that the structure suggested by the data of fig. 6 (as in fig. 3) is associated with 
Z-oscillations in dE/dx and in W[3S]. 

5. CONCLUSIONS 

Electronic sputtering and electron emission for keV ions incident on water ice are 
determined by the electronic stopping power and additional processes in the ionization track. The 
electronic sputtering yield grows linearly with the square of the electronic stopping power but the 
proportionality factor is velocity dependent for slow ions. The electron yield induced by light ions 
is larger for water ice than for metals, but the yield grows sub-linearly with electronic stopping 
power, due to the increase barrier produced by unbalanced holes in the insulator. A relatively 
simple model that includes a track potential acting as a spherical barrier can account for this 
effect. 
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FIGURE CAPTIONS 


Fig. 1. 
Fig. 2. 


Fig. 3. 


Fig. 4. 


Fig. 5. 


Fig. 6. 


UHV experimental apparatus. The polarity of the 300 V battery is 
reversed to measure sputter yields. 


Sputtering and electron yields of H 2 O ice films vs. proton energy. The 
hollow circle and filled triangle are 10-keV D + points. The squares are 
sputter yields from Brown et al. [20]. Only some representative relative 
error bars are shown, while the lines are to guide the eye. 


Electron yields and sputter yields from water ice as a function of the atomic 
number Z normalized to the proton yield; stopping power from Yarlagadda 
et al. [25] normalized to the proton stopping power. 


Electron yields from Al and water ice vs. projectile atomic number (Z) at 5 
keV/amu. Circles are yields from H 2 O ice films (this work), triangles are 
yields from Al (extrapolated from Alonso et al. [26]). 


The model for unbalanced holes forming an effective potential in a track. 
Outgoing electrons (• ) are attracted by the unbalanced holes (©) , 
which act to produce an additional barrier. 


Electron yield as a function of the electronic stopping cross section. The 
line is calculated using the model of Sec. 3C and is normalized to 1 .9 
electrons/ion at 10 eV/(10 15 H2O/cm 2 ). The experimental data is from the 
present work (proton data from fig. 2 and other ion data from fig. 4). 
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